Swelling of particle-encapsulating random manifolds 
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We study the statistical mechanics of a closed random manifold of fixed area and fluctuating 
volume, encapsulating a fixed number of noninteracting particles. Scaling analysis yields a unified 
description of such swollen manifolds, according to which the mean volume gradually increases 
with particle number, following a single scaling law. This is markedly different from the swelling 
under fixed pressure difference, where certain models exhibit criticality. We thereby indicate when 
the swelling due to encapsulated particles is thermodynamically inequivalent to that caused by 
fixed pressure. The general predictions are supported by Monte Carlo simulations of two particle- 
encapsulating model systems — a two-dimensional self-avoiding ring and a three-dimensional self- 
avoiding fluid vesicle. In the former the particle-induced swelling is thermodynamically equivalent 
to the pressure-induced one whereas in the latter it is not. 

PACS numbers: 87.16.D-, 64.60.Cn, 64.60.De, 68.35.Md 



I. INTRODUCTION 

There has been considerable interest in the past few 
decades in the statistical mechanics of membranes and 
surfaces [l[ . This has been partly motivated by the ubiq- 
uity of bilayer membrane vesicles [2J in various natural 
and industrial systems. Since the lateral size L of such 
envelopes is much larger than their thickness, they can be 
treated to a good approximation as (d — l)-dimensional 
manifolds, d being the embedding dimension. Another 
consequence of the thinness of the membrane is that it 
resists stretching much stronger than bending. Hence, 
the surface area A of the membrane is usually assumed 
fixed. The statistical mechanics of such a manifold in- 
volves an interplay between conformational fluctuations 
and bending elasticity, leading to a characteristic per- 
sistence length, l p [H, 0] ; over distances smaller than l p 
the manifold is essentially smooth, whereas beyond it the 
surface becomes random. When the manifold is closed (a 
vesicle), its smoothness is affected not only by the elastic 
persistence length but also by the degree of swelling (e.g., 
volume-to-area ratio). 

The various studies of vesicle thermodynamics can be 
classified in two groups according to the volume con- 
straint that they impose (for a given A). One body of 
works, e.g., Refs. [1, H, 0, Hi , considers the ensemble of 
fixed volume V. These studies, aimed at actual bilayer 
vesicles, assume the low-temperature limit, l p > L, in 
which the vesicle is represented by a continuous closed 
surface in three dimensions (3D). The various equilibrium 
shapes are derived as ground states of the elastic Helfrich 
Hamiltonian [tJ , depending on the dimensionless volume- 
to-area ratio, V/A z / 2 . Another body of works treats the 
ensemble of fixed pressure difference p across the man- 
ifold. The studied systems include Gaussian [§], freely 
jointed and self-avoiding [UdO! 

rings in two 
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dimensions (2D), as well as model fluid vesicles in 3D 
[H [H, [H M, ES HI IH- Most of these works assume 
the random, high-temperature limit (l p <C L), yet the 
crossover to l p > L was addressed as well [HI G~l> EH, ■ 

As long as equilibrium averages are concerned, the en- 
sembles of fixed V and fixed p are equivalent, i.e., they 
are related by a smooth, single-valued (Legendre) trans- 
form. We focus here on another swelling scenario, where 
the manifold encapsulates a fixed number Q of parti- 
cles while its volume is unconstrained. The interest in 
such particle-swollen manifolds is not merely theoretical; 
most actual vesicles are immersed in solution and their 
membrane, over sufficiently long, experimentally relevant 
time, is semipermeable, allowing solvent exchange while 
keeping the solute trapped inside [H, [23|, [24| • Note that 
the particle number Q does not a priori imply a certain 
osmotic pressure, because the manifold is free to change 
its mean volume and, hence, the mean particle concen- 
tration. Nonetheless, since the mean volume and pres- 
sure should monotonously increase with Q, one expects 
to find equivalence (i.e., certain well-behaved transforms) 
between the fixed-Q ensemble and the other two. We 
have recently demonstrated, however, that these ensem- 
bles are not equivalent for a freely jointed ring in 2D [25| . 
A second-order transition between crumpled and swollen 
states, which occurs in the fixed-p ensemble at a criti- 
cal pressure p c [l(|, disappears in the fixed-Q case. The 
criticality is avoided as the system selects such a mean 
volume that the mean pressure always lies above p c for 
any value of Q. Thus, the regions of phase space covered 
in the two ensembles are different. 

In the current work we generalize these results, obtain- 
ing a unified description for a (d— I)-dimensional random 
manifold in d dimensions, swollen by either a fixed pres- 
sure difference or a fixed number of trapped particles. We 
thereby clarify when the swelling scenarios are thermo- 
dynamically equivalent and when they are not. It should 
be borne in mind that, while the current work is focused 
on strongly fluctuating, random manifolds (L ^> l p ), the 
bending rigidity and size of real vesicles place them in 
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the low-temperature, smooth regime (L < l p ) (26|. 

We begin in Sec. |TT] with a heuristic scaling analysis, 
which nevertheless yields the correct qualitative swelling 
behavior as found in previously studied models. We then 
proceed to verify these general results using Monte Carlo 
(MC) simulations of two model manifolds: a self-avoiding 
ring in 2D (Sec. for which the two scenarios are 

found to be equivalent, and a self-avoiding fluid vesicle 
in 3D (Sec. IIV|> . for which they are not. The results are 
analyzed and summarized in Sec. [Vj 



II. SCALING ANALYSIS 

We apply a scaling theory [H,|23,[ll| to a closed (d-1)- 
dimensional random manifold, composed of N nodes and 
embedded in d dimensions. In response to perturbation 
(pressure difference p or Q noninteracting trapped parti- 
cles), the manifold is assumed to be divided into subunits, 
or blobs, containing g nodes each. The blobs are defined 
such that each of them stores a tensile energy equal to 
the thermal energy k^T = 1 [28[ , 

i^- 1 ~ 1, (i) 

where 7 is the surface tension induced in the manifold 
due to the perturbation, and is the projected area 
of a blob. At length scales smaller than the blob size £ the 
manifold is unaffected by the perturbation and assumed 
to obey the power law, 

Z*- 1 ~ 9", (2) 

where v is a swelling exponent characterizing the unper- 
turbed manifold statistics. At distances larger than £ the 
perturbation stretches the manifold. The total projected 
area is given by the number of blobs times the projected 
area per blob, 

R d - X ~ (N/g)^- 1 . (3) 

So far Eqs. (JTJ)— (J3J) have been independent of the nature 
of perturbation (p or Q). The difference between the 
two cases enters via the Laplace law, which takes the 
following forms in the fixed-p and fixed-Q ensembles, re- 
spectively: 

l/R ~ P, (4a) 
1 /R ~ Q/R d . (4b) 

Solution of Eqs. (fTj)-(|4ap leads to the following power 
laws for the fixed-p case: 

(V) ~ R d ~ N^(pN^)^=^, 

1 v(d-l) 

7 ~ (pN d - 1 ) d "- 1 . (5a) 

(This result, in a different form, has been already ob- 
tained in Ref. [HI].) Two observations readily follow 



from Eq. (|5"a[) . First, the characteristic pressure differ- 
ence, required to appreciably swell the manifold (i.e., to 
obtain R ~ N 1 ^ (d ~ 1 '>), scales as p ~ Af-VO*-!^ regard- 
less of v. This characteristic value reflects the interplay 
between the mechanical work of inflating an object of vol- 
ume ~ N d /( d ~ 1 \ and the surface entropy of ./V degrees of 
freedom, pN d ^ d ~^ ~ N. Second, in cases where dv = 1 
the exponents diverge, i .e.. the analysis breaks down, and 
one expects criticality [131 ]. Both conclusions are borne 
out by previously studied models. Gaussian 0] and freely 
jointed 10] rings, having d — 2 and v = 1/2 (i.e., dv = 1) 
behave critically at p c ~ N" 1 , the former swelling to in- 
finite volume, and the latter undergoing a second-order 
transition to a smooth state. By contrast, self-avoiding 
rings, w ith d = 2 and v — 3/4, swell gradually with p 

EMEUS- 

Turning to the fixed-Q case, we find from Eqs. (TTj)- (|4b[) 
the power laws, 

(V) ~ N^(Q/N) ± ^ 1 , 
7~(QW- (5b) 

The corresponding observations in this case are as fol- 
lows. First, appreciable swelling occurs for Q ~ N, 
regardless of v and d. Thus, the number of encapsu- 
lated particles required to swell the envelope scales with 
the area only, rather than the volume. This is a conse- 
quence of considering a vanishing external pressure [26j |. 
In such a case the particle entropy (~ Q) has to compete 
only with the surface one (~ N). Second, there is no 
divergence of exponents in Eq. (|5b[). i.e., no criticality. 
Both conclusions are consistent with findings regarding 
particle-encapsulating freely jointed rings in 2D [25j . 

The two blob analyses, along with the resulting power 
laws [Eqs. (f5aj) and (|5b[) ]. should hold so long as 1 < g < 
N. This corresponds to the restrictions, N~ dv ^ d ^ x ^ < 
p < N~ 1 \ and 1 < Q < N. At larger swelling, 
nonetheless, we expect the manifold to be smooth, having 
(V) ~ N d ^ d ~ 1 \ According to Laplace's law this leads to 
a surface tension 7 ~ pN 1 /^- 1 ) and 7 ~ Q/N. Combin- 
ing these large-swelling results with Eqs. l[5"a]) and {5b}, 
and provided there is no criticality (dv ^ 1), we conjec- 
ture the following scaling relations, expected to hold for 
all values of p and Q : 

(V) = N^f p (pN^), 

7 = h p ( P N^), (6a) 

(V) = N^f Q (Q/N), 
7 = h Q (Q/N). (6b) 

The scaling functions for the mean volume, f p and fg, 
should cross over from the power laws of Eqs. (|5a|) and 
(|5b[) for small arguments to constant values for large ar- 
guments. The scaling functions for the surface tension, 
h p and Hq, are expected to cross over from the power 
laws of Eqs. (|5a|) and (|5b[) to linear ones. The valid- 
ity of Eqs. (|5bp and (|6b[) has been already proven for a 
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particle-encapsulating freely jointed ring in 2D [25|. In 
addition, the scaling of Eq. (jrJa|) has been demonstrated 
in the swelling of those rings with increasing p above the 
critical point [H|]. We now proceed to check the validity 
of Eqs. (5)-(6) in two additional model systems. 



III. SELF-AVOIDING RING IN 2D 

We first follow the model and MC scheme presented in 
Refs. [H, [TH, [3] for a 2D self-avoiding ring subject to 
an inflating pressure difference p . The manifold is repre- 
sented by a closed chain of N self-avoiding circles (beads) 
of diameter a = (5/9)1, linked by tethers of maximum 
length 1 = 1. In each MC step every bead is moved to a 
random position within a square of (—0.2(7, 0.2<t) 2 about 
its former position. These values of u, I, and maximum 
step size ensure that self-intersection of the ring cannot 
occur. The move is weighted by W = e pAV , where A.V 
is the difference in (2D) volume of the ring due to the 
move, and is accepted provided that (i) self-avoidance is 
fulfilled; (ii) tethers do not exceed their maximum length; 
and (iii) W exceeds a random number in the range [0,1]. 
Simulations were performed for TV = 50 to 800. 

The mean volume of the ring as a function of pres- 
sure difference is presented in Fig. [1] The different data 
sets collapse onto a single curve once the mean vol- 
ume is rescaled by the maximum volume of the ring, 
Vmax = N 2 /(A-k), and the pressure by iV -1 , in accord 
with Eq. (|6a|) . This scaling law, however, yields a van- 
ishing mean volume for p = 0, whereas the unperturbed 
ring has a finite mean volume of Vq ~ iV 2l/ , v = 3/4. 
In the thermodynamic limit (N — > oo) the correction is 
negligible, Vb/Vmax ~ N~ x /' 2 — > 0, but for finite rings the 
scaling of Eq. ((6aj) breaks down for sufficiently small p, 
as seen in Fig. [1] Therefore, to capture the initial lin- 
ear dependence of (V) on p, as predicted by Eq. (jS"a|) for 
d = 2 and v = 3/4, we replot in Fig. [T] (inset) the data for 
(V) — Vq. The initial increase of (V) — Vq with p seems 
to be consistent with a linear law, although we cannot 
claim to have clearly confirmed it. 

Next, we turn to particle-encapsulating manifolds by 
setting p — and introducing Q ideal particles at ran- 
dom positions inside the ring. Hard-core repulsion is in- 
troduced between the particles and envelope beads (but 
not between the particles themselves) , with particle-bead 
minimum distance of a. The MC step is extended to 
include repositioning of each particle within a square 
of (—0.2(7, 0.2er) 2 about its former position. This maxi- 
mum step size, together with the hard-core repulsion be- 
tween particles and envelope beads and maximum tether 
length, ensure that particles cannot exit the ring. Rings 
of N = 50 to 800 have been simulated, with Q ranging 
between and 20iV. 

In Fig. [2] we present the mean volume as a function 
of Q. In agreement with Eq. (|6b[) all data collapse onto 
a single curve when (V) is scaled by TV 2 and Q by N. 
As in the case of fixed p, discussed above, scaling breaks 
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FIG. 1: (Color online) Mean volume of 2D self-avoiding rings 
as a function of pressure difference. Data were obtained by 
MC simulations for different ring sizes N and rescaled accord- 
ing to Eq. (|6ap . V ma x = N 2 /(An) being the maximum volume 
of the ring. Inset shows the same data on a log-log scale after 
the mean volume of the unperturbed ring, Vq ~ N 3 ^ 2 , has 
been subtracted from (V). A solid line of slope 1 is shown for 
reference. 



down for very small Q, when (V) becomes affected by the 
finite volume of the unperturbed state. The power law 
predicted by Eq. (|5b[) for 2D self-avoiding rings (d = 2, 
v = 3/4), (V) ~ V^x (<2 /AO 1/2 , is nevertheless verified 
after subtracting Vb from the mean volume (Fig.[5]insct). 

To demonstrate the equivalence of the fixed-p and 
fixed-Q scenarios for this system we transform the data 
for pressurized rings (Fig. [I]) according to Q(p) = p(V(p)) 
and present them in Fig. [J] alongside the data for fixed 
Q. The data sets of the two scenarios match nicely over 
the entire ranges of p and Q. 



IV. SELF-AVOIDING FLUID VESICLE IN 3D 

The second manifold we consider is a discrete model 
of a fluid vesicle, which was extensively studied by MC 
simulations under fixed pressure difference p [H, [3, U3j 
[l8l . Il9l [2pj . The vesicle is represented by a closed, tri- 
angulated, off-lattice network of iV nodes (self-avoiding 
spheres) of diameter a = l/y2, interconnected by a fixed 
number of tethers of maximum length 1 = 1. Membrane 
fluidity is mimicked by constantly varying the network 
connectivity. The MC step comprises two parts, (i) Each 
bead is moved randomly within a cube of (— O.ler, 0.1c) 3 
about its former position (self-avoidance permitting). 
The move is weighted by a Boltzmann factor of e pAV , 
where AV is the change in volume caused by the move, 
(ii) iV attempts are made to break a randomly chosen 
tether, which has formed the common side of two tri- 
angles, and rebuild it between the two other corners of 
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FIG. 2: (Color online) Mean volume of 2D self-avoiding rings 
as a function of number of encapsulated particles. Data were 
obtained by MC simulations for different ring sizes N and 
rescaled according to Eq. ()6b[l . Vm ax = N 2 /(4ir) being the 
maximum volume of the ring. Also plotted are the data points 
from the fixed-p simulation (Fig.[l|, whose horizontal coordi- 
nate is calculated as p{V(p)) /N . Inset shows the data on a 
log-log scale after the mean volume of the unperturbed ring, 
Vo ~ N 3 / 2 , has been subtracted from (V). A solid line of 
slope 1/2 is shown for reference. 



those triangles (provided that the required tether length 
does not exceed I). The choice of a, I, and maximum 
step size prevents a bead from passing through another 
part of the network, making the manifold self-avoiding. 

The swelling of this model vesicle as a function of p fol- 
lows three regimes [lj| . (i) At low pressures the vesicle is 
in a collapsed state, having branched-polymer statistics, 
where the mean volume and mean-square radius of gyra- 
tion scale as (V) ~ R 2 ~ N, with negligible dependence 
on p [TH, HH, Ho, El . (ii) At a critical pressure, p — p* (N) , 
the vesicle undergoes a first-order transition to a swollen 
state, whose mean volume gradually increases with p as 
(V) ~ pO.47jyi.73 [15^ ( m ) ^ t sufficiently large p the 

power-law behavior crosses over to asymptotic swelling 
toward the maximum volume. 

The blob analysis presented in Sec. [Til obviously fails 
in regime (i) of low swelling, since the volume enclosed 
in such collapsed manifolds does not follow the standard 
relation (V) ~ R d . Instead, we use the fact that the ratio 
between the cross-section (frame) area of the manifold 
and its real surface area is vanishingly small. Such a 
manifold has a constant surface tension, 7 *~ 1 (in units of 
k B T/l 2 Applying Laplace's law, p ~ 7/i? ~ N^ 1 / 2 , 
we find that the deflated regime (i) is valid for p < N~ x / 2 , 
i.e., 
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(V) ~V ~N. 



(J) 



In regime (ii) the scaling analysis of Sec. HT1 holds. Com- 
parison of the previously obtained power law, (V) ~ 



f-^N 1 - 73 , with Eq. (HaJ for d = 3, gives v = 0.787 
fla ]. Our modified scaling analysis [Eqs. (f5"aj) and (jfja)) ] 
indicates that p is scaled with TV^ 1 / 2 . We note that the 
power-law dependence of the critical pressure p* on N 
has been controversial p^ |. with exponents ranging be- 
tween —0.5 [l|| and —0.69 [2(j. The scaling argument 
for p > p*, together with Eq. ([7]) for p < p* , strongly 
suggest that p* ~ N^ 1 / 2 [33[. 

We have repeated the MC simulations for fixed p, as 
presented in Refs. [lj|, while extending them to larger 
vesicles and higher pressure values. The results are shown 
in Fig. O scaled according to Eq. (|6"a|) . The first-order 
transition at p* ~ iV -1 / 2 is clearly reproduced, and the 
predicted scaling for the entire range of p > p* is con- 
firmed. The scaling for p > p* is not inconsistent with 



the power law of Ref. [15| and Eq. (|5a[) , having v between 
0.7 and 0.8 (Fig. [3] inset), yet this regime is too narrow 
to be clearly resolved. 
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FIG. 3: (Color online) Mean volume of 3D self-avoiding fluid 
vesicles as a function of pressure difference, as obtained by 
MC simulations for different vesicle sizes N. Data is scaled 
according to Eq. (I6a|) . exhibiting a discontinuous transition 
at p* ~ 7V _1//2 . For p > p* the data collapse onto a single 
curve. Inset presents the same data for p > p* on a log-log 
scale. A solid line of slope 0.6 (corresponding to v = 0.75) is 
shown for reference. 

We now turn to particle-encapsulating vesicles. Re- 
peating the aforementioned argument for the deflated 
state of regime (i), Q/(V) ~ j/R ~ iV -1 / 2 , we find 



(V) ~ N 3 ' 2 {Q/N). 



(8) 



(This linear dependence of (V) on Q will be shown in 
Sec. [V] to be intimately related to the phase transition 
observed as a function of p.) The scaling law of Eq. © 
for the low-swelling regime turns out to be similar to that 
of Eq. (|5b[) . Hence, despite the inadequacy of the blob 
analysis in regime (i), we expect the scaling conjecture, 
Eq. (|6b[) , to hold for all values of Q in this model as well. 
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To check these predictions we modified the MC scheme 
presented above by setting p — and adding Q ideal 
particles of diameter a, randomly positioned inside the 
vesicle. The particles do not interact with each other 
but have a hard-core repulsion with the network nodes, 
keeping them trapped inside the vesicle. The MC step is 
extended to include random repositioning of each particle 
within a cube of (— 0.1c, O.lcr) 3 about its former position. 
Vesicles with N ranging between 162 and 642 and Q up 
to ION (for the smallest vesicle) have been simulated. 

Results for the mean volume as a function of Q for 
various vesicle sizes are shown in Fig. [4] Once the volume 
Vb of the unperturbed (branched) state [Eq. (J7J)], which 
is inaccessible to particles due to their excluded-volume 
interaction with the manifold, is subtracted from (V), the 
data collapse onto a single curve according to Eq. (|6b[) . 
Two power-law regimes are seen in Fig. 0] (inset). At 
low swelling (V) increases linearly with Q, in agreement 
with Eq. © 0. At about Q ~ 0.087V the swelling 
crosses over to a different power law which, when fitted 
to Eq. ([5b]). yields v = 0.75(2). This value is close to that 
found in the fixed-p simulations, v = 0.787 [l5ll35j. For 
larger values of Q this power-law regime should cross over 
to asymptotic saturation toward the maximum volume. 
Because of computer limitations we could sample only 
the lowest edge of this regime (Fig. [4] inset). 



0.04 



' * 


N = 


162 


X 


N = 


252 


- o 


N = 


362 


" □ 


N = 


492 


o 


N = 


642 











0.001 t 



12 3 4 

Q/N 

FIG. 4: (Color online) Mean volume of 3D self-avoiding 
fluid vesicles as a function of the number of trapped par- 
ticles, as obtained by MC simulations for different vesicle 
sizes N. Data collapse onto a single curve according to Eq. 
(|6b|) once the volume of the unperturbed vesicle, Vb ~ N, 
is subtracted from (V). Inset represents the same data 
on a log-log scale, exhibiting a linear regime for Q <C N , 
((V) - V )/N 3/2 ~ (Q/iV) 102(3) (dotted line), followed by a 
more swollen regime with {(V) - V )/N 3/2 ~ (Q/iV) a38(3) 
(dashed line). The arrow indicates the crossover between the 
two regimes at Q ~ 0.08AT. 

Unlike the case of fixed p, the vesicle gradually swells 
with Q, exhibiting no phase transition. To further ver- 



ify the absence of a first-order transition we have mea- 
sured the probability distribution function of the vol- 
ume, P(V), as a function of Q. Whereas under fixed 
p, at p — p* , one finds a bimodal distribution [Ref. (l5| 
and Fig. |5(a)], i.e., coexistence of collapsed and swollen 
states, for particle-encapsulating; vesicles we obtain uni- 
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FIG. 5: (Color online) Volume distribution functions of 362- 
node 3D vesicles, (a) For fixed pressure difference, at the 
transition point, the distribution is bimodal. (b) Under fixed 
number of encapsulated particles the distribution is always 
unimodal. 

Finally, let us consider the effective pressure exerted 
by the encapsulated ideal particles, p — Q/(V). From 
Eq. (|6b[) we have 



p = N-^iQ/N), Hx) = x/f Q (x). 



(9) 



In the low-swelling regime we have found a linear behav- 
ior, f Q (x < 1) ~ x, [Eq. (0) and Fig. H inset]. Thus, 
ip(x < 1) = const, i.e., the effective pressure does not 
change with Q throughout this regime. Figure [5] demon- 
strates the data collapse according to Eq. ([5]), as well 



6 



as the finite, constant pressure p m in at low swelling even 
for the smallest values of Q. (In calculating the concen- 
tration and pressure from the simulations we have con- 
sidered the particle- accessible volume, (V) — Vb.) One 
expects Pmin to coincide with the transition value under 
fixed pressure, p*. [Compare also Fig. O^a), plotted for 
p = p* , with Fig. G3h) , where the effective pressure is es- 
sentially fixed at p m in for all curves.] We find, however, 
p* ~ 1.8p m in- This discrepancy may stem from the in- 
teraction of the particles with the vesicle, making them 
deviate from the ideal-gas behavior, particularly in the 
deflated state. 



18 . * N= 162 
^ - >< N=252 
k,°, 6 : o N=362 
. - □ /V = 492 
> ; o /V=642 
^ 14- 



O) 12- 
II 



> 10 

ft, 



o.i 



Q/N 



0.2 



0.3 



FIG. 6: (Color online) Effective pressure of encapsulated par- 
ticles as a function of particle number, as obtained by MC 
simulations for various vesicle sizes N. Data collapse accord- 
ing to Eq. |9j. Note the finite effective pressure at small Q. 
The arrow indicates the crossover between the two swelling 
regimes at Q ~ 0.087V. 



V. DISCUSSION 

The scaling analysis presented in Sec.|TT]yields a unified 
account of the swelling of random manifolds with increas- 
ing pressure difference or number of encapsulated parti- 
cles. The validity of this description has been demon- 
strated for several model systems in Sees. IIIII and IIVI 
and in Ref. [25(. Similar scaling analyses for the case of 
fixed p were previously presented in Refs. [H, [l5| . Those 
analyses and ours coincide in the power-law regime, Eq. 
(I5a[) . However, while the previous analyses are focused on 
the weak-swelling regime and constructed to include the 
random, unperturbed state of the manifold, the one pre- 
sented here is aimed at encompassing the high-swelling 
behavior. Thus, on the one hand, our scaling relations, 
Eqs. (|6aj) and (fob)) , cannot account for the unperturbed 
state and give a vanishing mean volume in the limit of 
vanishing perturbation. The range of p (or Q) where this 
deficiency is relevant, nonetheless, vanishes in the ther- 



modynamic limit [361 ]. On the other hand, whereas the 
previous analyses assumed that scaling broke down at 
sufficiently large swelling [ID, [TBI , we have claimed that 
Eqs. (|6a[) and (|6b[) should hold for the entire range of 
p or Q. Although there is a priori no reason why the 
scaling behavior should have this broad range, we an- 
alytically proved the conjecture for freely jointed rings 
in 2D at fixed p > p c [lCj] or fixed Q [25[ . In the cur- 
rent work we have provided further numerical support 
of the scaling conjecture in several additional systems — 
2D self-avoiding rings at fixed p or fixed Q, and 3D fluid 
vesicles at fixed p > p* or fixed Q. Hence, provided that 
the swelling exhibits no criticality, the scaling relations 
Eqs. (|oa)l and (fob)) seem to be applicable in a broad range 
of systems. 

The analysis presented here for the swelling with Q 
yields new scaling relations, which have been confirmed 
in all studied systems. The different behavior of particle- 
encapsulating manifolds lies in the response of manifold 
volume to changes in the number of encapsulated parti- 
cles. The resulting particle concentration and effective 
pressure depend on this response and, therefore, may 
have a nontrivial dependence on Q. In certain cases this 
may lead to thermodynamic inequivalence of the fixed-p 
and fixed-Q ensembles. Equivalence breaks down when 
the two ensembles are no longer related by a one-to-one 
smooth transform. In the two examples where inequiv- 
alence has been demonstrated — 2D freely jointed rings 
[25| and 3D fluid vesicles (Sec. IIV[) — both conditions 
of smoothness and single-valuedness are violated: (i) a 
criticality in the fixed-p ensemble makes the transform 
Q(p) = p(V(p)) nonanalytic; and (ii) the effective pres- 
sure in the fixed-Q ensemble is bounded from below by 
a finite value, i.e., states of low pressure are inaccessible 
(cf. Fig. [S]). We now show that this combination of criti- 
cality under fixed p and inaccessible states for fixed Q is 
not a coincidence. 

Let us consider a general power-law response to par- 
ticle number, (V) ~ Q a . Transforming to the fixed-p 
ensemble, we get p{Q) = Q/(V) ~ Q^ a and (V) ~ 
pa/(i-a)^ Several observations follow from these rela- 
tions. First, thermodynamic stability dictates that (V) 
increase with p, i.e., a < 1. We are left with two different 
cases, (i) If a < 1, there is no criticality and arbitrar- 
ily small values of Q will correspond to arbitrarily small 
values of p. Hence, in this case there is equivalence, (ii) 
If a = 1, we expect both criticality under fixed p and 
inaccessibility of small-pressure states at fixed Q, i.e., 
inequivalence of the two swelling scenarios. For maxi- 
mum a the manifold volume is maximally susceptible to 
changes in Q (linear in Q), to the extent that the con- 
centration and pressure do not change with Q (cf. Fig. 
[5]). Thus, criticality and inequivalent phase spaces come 
hand in hand. In a standard case, where the blob analysis 
of Sec.|H]holds, we get from Eq. j5b} a = d(l — u)/(d—l), 
and the condition a < 1 is equivalent to dv > 1. In addi- 
tion, one has a geometrical lower bound for the swelling 
exponent, which cannot be smaller than that of a folded, 
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compact manifold, v > (d—l)/d. This leads to the re- 
striction a < l/(d—l), which is consistent with, and 
stricter than, the thermodynamic one, a < 1. Hence, 
we conclude that for most systems, which obey the anal- 
ysis of Sec. im case (ii) above, involving criticality and 
inequivalence, can occur only in 2D, i.e., for d = 2 and 

i/= 1/2 Gsm. 

All of these general conclusions are supported by spe- 
cific examples. A 2D self-avoiding ring is an example of 
case (i) above. It obeys the scaling analysis of Sec. HTl with 
d = 2, v = 3/4, i.e., a — 1/2. (The value of a has been 
confirmed by simulations; see Fig. [51) This system ex- 
hibits no criticality under fixed p, and the two ensembles 
have been found equivalent (Fig. (2). The more interesting 
case (ii) has been encountered in three systems. Two ex- 
amples are provided by Gaussian and freely jointed rings 
in 2D [23|. For both examples the blob analysis holds, 
and d = 2,v = 1/2 [i.e., dv = 1). The third example 
of the anomalous case (ii) is a 3D fluid vesicle fSec. IIV|) . 
for which the blob analysis of Sec. [Til fails, yet a linear 
dependence of (V) on Q has been found [Eq. © and Fig. 
3]. Indeed, under fixed p 0, [T3, EH all three examples 
exhibit phase transitions. 

This host of examples leads to the expectation that 
the picture described here, including the scaling relations 
and possible phase transitions, should hold for any ran- 
dom manifold swollen by either a pressure difference or 
encapsulated particles. In cases where the blob analysis 
of Sec. [TIJis valid, one needs to know merely the dimen- 
sionality d and the statistics of the unperturbed manifold 
(v) to predict the qualitative swelling behavior. In other, 
exceptional cases (e.g., the 3D fluid vesicle of Sec. IIV[) it 
suffices to know the response of the unperturbed man- 
ifold to a small number of encapsulated particles (i.e., 
a). 

The thermodynamic inequivalence between the fixed- 
p and fixed-Q scenarios, reported above for certain sys- 
tems, also implies inequivalence between the canonical 
and grand-canonical ensembles in those systems. This is 
because fixing the chemical potential (i of the encapsu- 
lated particles inevitably fixes also the mean pressure p 
that they exert on the manifold, as these two intensive 



variables are related via the particles' equation of state. 
(For example, for ideal particles \i = lnp.) Once again, 
because of the unconstrained volume, the system cannot 
attain arbitrarily small concentrations as Q is decreased, 
and, therefore, the full range of /i is not covered. The 
inequivalence of the fixed-Q and fixed-// ensembles has 
been directly demonstrated for freely jointed 2D rings 




In the current work we have not explicitly considered 
the bending rigidity of the vesicle. Such a bending-free 
description is valid in two limits: (i) At sufficiently strong 
swelling the fluctuations of any vesicle are governed by 
surface tension rather than bending rigidity, (ii) If the 
manifold is sufficiently large (L 3> £ p ), bending rigidity 
merely renormalizes the molecular length a to l p and the 
number of surface degrees of freedom N to Na 2 /lp. It 
is this random, strongly fluctuating case which has been 
the focus of the current work. On the one hand, due to 
their bending rigidity (k ~ 10 /cb? 1 ) and size (0.1-10 /im), 
real bilayer vesicles are smooth and do not satisfy limit 
(ii). On the other hand, some of our most significant re- 
sults (e.g., the inaccessibility of low-pressure states) con- 
cern weak swelling, outside limit (i), where the bending 
rigidity of real vesicles plays an important role. Thus, 
the direct relevance of the current work to real bilayer 
vesicles is limited. Yet, overall, this work and the spe- 
cific examples associated with it highlight the qualitative 
differences which may emerge between pressurized mani- 
folds and particle-encapsulating ones. Indeed, the differ- 
ent behavior of particle-encapsulating vesicles is manifest 
also in realistic scenarios involving smooth membranes, 
e.g., highly swollen bilayer vesicles in solution [26| . 
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